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The doubly heavy baryon Hqq/ (Q''' = ti or c) is different from the ordinary baryons. The 
production of the doubly heavy baryon can provide valuable insight on how the colored {QQ')- 
diquark can be transformed into the color-singlet baryon. We study the doubly heavy baryon 
production through the e"*"e~ annihilation channel, e"*" -|- e~ — >■ 'y/Z° — >■ Sqqi + Q + Q', within the 
nonrelativistic QCD framework. We calculate the baryon transverse momentum and the rapidity 
distributions for all these channels. Typical transverse momentum and rapidity cuts are adopted 
to show the properties of these distributions in detail. At a e'^e~ collider that runs around the 
Z^-boson mass with a high luminosity up to £ ~ 10^*~^®cm~^s~^, in comparison to the Belle and 
BABAR experiments, it is found that sizable Hcc, Sbc and Hi,;, events can be produced even after 
performing reasonable baryon transverse momentum and rapidity cuts. 
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Within the nonrelativistic QCD (NRQCD) frame- 
work the physical state of a heavy quarkonium is 
generally regarded as a superposition of Fock states, and 
the relative importance among those infinite ingredients 
is evaluated by the velocity scaling rule. The production 
of the doubly heavy baryon Sqq/ can provide valuable 
insight on how the colored (Q(5')-diquark can be trans- 
formed into the color-singlet baryon. These baryons offer 
a good platform for testing various theories and models. 
Throughout the paper, if not specially stated, we shall 
ignore the isospin-breaking effect of the doubly heavy 
baryons; namely, the symbol Sqq' is a short notation for 
all the baryon Sqq/^, where Q^'-* — b or c quark and the 
light quark q stands for u, d or s, respectively. 

In Ref . 0] we have calculated the doubly heavy baryon 
production through the e'^e~ annihilation, + e~ 
7/ZO Eqq,+Q+Q\ within the framework of NRQCD. 
We have discussed the total cross section for the doubly 
heavy baryon production through e"'"e~ annihilation at 
the super-Z factory a e+e^ colhder with a high lu- 
minosity C (X 10'^*~'^^cm~^s~^ and its colliding energy 
is around the peak. Moreover, the International Lin- 
ear Collider (ILC) is also programmed to run at the Z° 
peak with a luminosity of £ = 0.7 x 10"^*cm~^s~^ (GigaZ 
program [J). It has been observed that the Z°-boson 
resonance effect shall raise the production rate by sev- 
eral orders of magnitude in comparison to the Belle and 
BABAR experiments, then it could be a better platform 
for studying the properties of heavy baryons. It is in- 
teresting to show more detailed production properties, 
such as how the baryon production cross section varies 
with its transverse momentum or rapidity, which is ex- 
perimentally more useful. 



According to NRQCD factorization theorem [H, [H, the 
differential cross section of the e+e^ annihilation process 
takes the following factorization form. 



dcr(e^ 



(QQ')[n] + Q- 



Q') 

Q'){0"{n)), (1) 



where the matrix element {O^ (n)) is proportional to 
the inclusive transition probability of the perturbative 
{QQ')[n\ pair into the doubly heavy baryon Sqq/, and 
the symbol [n] represents the spin and color quantum 
numbers for the diquark state (QQ'). Since the fragmen- 
tation function Djz) for a heavy diquark into a baryon 
peaks near z « 1 [6|-l8| , the momentum of the final baryon 
may be considered roughly equal to the momentum of 
initial diquark. Thus, as an approximation, studying the 
hadronic production of Sgg/ is equivalent to studying 
the hadronic production of (Q(5')-diquark. However, dif- 
ferent to the heavy quark fragmentation, during the frag- 
mentation of a diquark into a baryon, the diquark may 
dissociate Q, which will decrease the baryon production 
cross section to a certain degree. In the present paper, 
we will not take such dissociation effect into considera- 
tion. Regarding this point, our present estimations can 
be treated as an upper limit for the production chan- 
nels. For the intermediate (cc)- and (bb)- diquarks, there 
are two spin and color configurations: pS'ija and [^S'oje; 
while for the (6c)-diquark, there are four spin and color 
configurations: [^Si]^, [^Si]e, [^^olsi and [^5o]6- 



* wuxg@cqu.ed u. cii] 
t |smwang@cqu.edu.cn| 



^ Due to the heavy quark-diquark symmetry I9l 4llll . there are re- 
lations between the doubly heavy baryons and the singly heavy 
mesons, e.g. they could have similar properties such as the spcc- 
trums and the productions or the decays. Then the doubly heavy 
diquarks may dissociate due to the radiation of soft partons with 
large probability as those of the heavy mesons. 
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TABLE I. Total cross sections (in unit: pb) for the HnQ'(n) 
baryon production through e^e^ annihilation at the Z peak 
{y/S = mz) or the Higgs peak (v^ = m_r/) respectively, where 
n stands for the intermediate (Q(5')-diquark state. 





<j{y/S — mz) 


<j{VS — mn) 


e+e 7 Scc(pS'i]3) 


8.90 X 10"'' 


4.81 X 10-" 


e+e" -)■ 7 Scc([^So]6) 


4.29 X 10-'* 


2.32 X 10-" 


e+e- ^ ^« -> Hec(pSi]3) 


0.44 


8.12 X 10-" 


e+e- Z° ^ Hccd'Sole) 


0.21 


3.94 X 10-" 


e+e- 7 H6c(p5'i]3) 


2.66 X 10-" 


1.45 X 10-" 


e+e- 7 Hbc(p5i]6) 


1.33 X 10-'* 


7.24 X 10-^ 


e+e- 7 H6c([^5o]3) 


1.85 X 10-" 


1.05 X 10-" 


e+e- -s> 7 -s> Sbc([^5o]6) 


9.27 X 10-^ 


5.27 X 10-^ 


e+e- ^ Z« ^ H,,(pSil3) 


0.61 


1.15 X 10-^ 


e+e- ^ ^0 ^ EU['Si]e) 


0.30 


5.73 X 10-" 


e+e- ^ ^ Si,c(rSo]3) 


0.44 


8.43 X 10-" 


e+e- ^ H6c(['5ol6) 


0.22 


4.21 X 10-" 


e+e- 7 ^ H6b(p5'i]3) 


1.94 X 10-^ 


1.12 X 10-^ 


e+e- 7 H(,5([^5o]6) 


8.95 X 10-" 


5.28 X 10-*^ 


e+e- ^Z°^ Ebt{('Si]s) 


4.85 X 10-^ 


9.62 X 10-^ 


e+e- ^ Z" ^ Hi,i,([iS'o]6) 


2.33 X 10-^ 


4.64 X 10-^ 



The short-distance cross section 
da (e+e- ^ (QQ')W + Q + Q') = P^^'f^' - 

where M is the hard scattering amplitude, and ^ means 
we need to average over the spin states of the electron 
and positron and the sum over the color and spin of all 
final particles. The detailed treatment of M has been 
presented in Ref. Q , the key point of which is to properly 
transform one of the fcrmion lines of the diquark produc- 
tion to an antifermion line of the meson production. d$3 
stands for the conventional three-particle phase space, 
which can be generated and integrated totally or partly 
by using the FormCalc program [l2| , or a combination of 
RAMBOS [11 and VEGASjll, which can be found in 
the generators GENXICC [ijand BCVEGPY 

When doing the numerical calculation, we take the in- 
put parameters as = 2.4952 GeV, mz = 91.1876 
GeV, mw = 80.385 GeV and mu = 125 GeV. The renor- 
malization scale for Sqq' as well as other input param- 
eters, such as the bound state parameters and the non- 
perturbativc matrix elements, arc taken to be exactly the 
same as those adopted in Ref. [3]. 

Total cross sections for the Sqq' in) baryon production 
through the e+e- annihilation at the peak {\fS = 
mz) or the Higgs peak {\/S = mn) are presented in 
Table H] respectively, where n stands for the intermediate 
((5(3')-diquark state 0. 



^ We correct an input error for the parameter g„ used in Ref. 
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Zo,n = PSi]5 




Zo,n = i'Soja 











10' 10^ 10^ 



E GeV 

cm 

FIG. 1. Total cross section (in unit: pb) for the Hcc- 
production channel, e+ 4- e- — >■ 7/^" — >■ 'E.cc{n) -I- c-(- c, ver- 
sus the e+e- collision energy Ecm = "s/S, where n stands 
for the intermediate (cc)-diquark state. Here the peak is at 
y/S — mz- 

When the e+e" collider runs at the peak {^/S = 
mz), the total cross section a{mz) for the production 
channel via the 7-propagator is about three orders lower 
than that of the channel via the Z°-propagator. For the 
collider runs at the Higgs peak {^/S — mn), the total 
cross section aimn) via the 7-propagator is suppressed 
by only one order compared to that of the channel via the 
Z'^-propagator. As an explicit example, we present the 
total cross section for the production channel e+ -I- e- — >■ 
j/Z'^ — > 'E.cci'n) + c + c versus the e+e- collision energy 
Ecm = in Fig.dl]). Moreover, we have 

• By summing up all the intermediate (Q(5')-diquark 
states' contributions into consideration, we obtain 
(Jtotimz) = 0.65 pb and atotimR) = 1.92 x 10-^ 
pb for Sec production; atot{mz) = 1.57 pb and 
o'totinT'H) = 3.36 X IQ-^ pb for S^c production; and 
o-tot{mz) = 0.07 pb and (Jtot{mH) = 1-59 x lO"'' 
pb for S^b production. 

• In contrast to the BABAR and Belle, at the super- 
Z factory, the production cross section for S^c is 
about two times larger than that of Sec- Thus, the 
super Z factory shall be helpful for studying the 
properties of Sbc- 

• For the production channels through the 7- 
propagator, the cross sections at the Z^ peak and 
the Higgs peak are on the same order of magnitude, 
i.e., about lO-^-lO-"' pb, which are small and can 
be safely neglected in comparison to the production 
channel through the Z'^-propagator. 

• As for the production through the Z'^-propagator, 
the cross sections a^mn) are almost three orders 



which should be gw = 2\/2mw y Gp/\/2- Then, all the numer- 
ical results listed in Ref. Q| for the total cross-sections via the 
Z'^-propagator should be corrected and multiplied by an overall 
parameter (mw/m^)* = 0.6039. All the Fortran sources are 
available upon request. 
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smaller than a(mz) for all the mentioned doubly 
heavy baryons. This shows the super-Z factory is 
much better than the super-Higgs factory for study- 
ing the doubly heavy baryon properties. 



TABLE II. The events of Hqq/ (n) through the production 
channel e^e~ — >■ 'y/Z'^ — >■ Hqq/ (n) + Q + Q' at various exper- 
iments, where n stands for the intermediate diquark state. 





super-.Z 


GigaZ 


supcr-Higgs 


Hcc(PS'l]3) 


4.40 X 10" 


3.08 X 10"* 


129.3 


Hcc(["''S'o]6) 


2.14 X lO"* 


1.50 X 10" 


62.6 


Hbc(PSl]3) 


6.07 X 10'' 


4.25 X 10" 


129.1 


HfcclP Si]e) 


3.03 X lO"* 


2.12 X 10" 


64.5 


S.bc{[^So]3) 


4.41 X 10'' 


3.09 X 10" 


94.8 


Ebcil^ So]e) 


2.21 X lO"* 


1.55 X 10" 


47.4 


H..(pSi]3) 


4.86 X 10^ 


3.40 X 10^ 


10.7 


S.bbil^ So]e) 


2.33 X 10^ 


1.63 X 10^ 


5.2 



Considering the super-Z factory and the GigaZ pro- 
gram of ILC, we can estimate the Sqq' baryon events 
generated per year. Here, we adopt the luminosity of 
super-Z factory as £ (x 10'^''^cm~^s~^. Our estimations 
are put in Table |TT1 where as a comparison, we also 
present the events at the super-Higgs factory 0. Ta- 
ble Hi] shows that both the super Z factory and the GigaZ 
program can generate considerable baryon events, e.g. 
^ 6.5 X 10^ Sec events, ~ 1.6 x 10^ S^c events and 
~ 7.2 X 10'^ events can be generated at the super 
Z factory; and the events per year under GigaZ are 0.7 
times of that of the super-Z factory. Such a large number 
of events shall provide us great chances to study the prop- 
erties of the doubly heavy baryons, such as to discuss its 
distributions and to further distinguish the baryons with 
different light quark contents, i.e., to distinguish Sqq'„, 
^QQ'd and rigg's- 

Next, we discuss the baryon transverse momentum (pt) 
and rapidity (y) distributions for \/S = mz- We present 
the pt- and y- distributions for the production channels 
e+ + e- ^ 7/ZO ^ Eqq, (n) + Q + Q' in Figs. (l2l3l4l5L 
which are for ^ccin), ^bb{n) and ^bc{n) respectively. 

Experimentally, there is no detector that covers all the 
kinematics, only some of the events can be measured. 
Especially, the baryon events with a small pt and/or a 
large rapidity y, i.e. the produced baryons moving close 
to the beam direction, cannot be detected directly, so 
such kind of events cannot be utilized for experimental 
studies. Therefore, the baryon events with proper kine- 
matic cuts on Pt and y must be put on precisely in the 



^ a e+e~ collider running at the Higgs peak with the same lu- 
minosity as that of the super-Z factory, which has recently been 
suggested and stimulated by the discovery of the standard model 
Higgs-like particle at the LHC. 




-3-2-10123 



y 

FIG. 2. The pt~ and j/— distributions of the production 
channel e+ 4- — >■ 7/^" — >■ 'S.cc{n) -f c -f c at peak 
= mz), where n stands for the corresponding interme- 
diate (cc)-diquark state. 



estimates. Considering detectors' abilities and to offer 
experimental references, we try various cuts in estimat- 
ing the baryon production. We introduce the parameter 
Ptcut with ptcut = A means pt > A and the parameter 
ycut with ycut = B means \y\< B. 



TABLE III. Total cross sections (in unit: pb) at the Z'^ peak 
{■\fS — mz) for the Hqq/ baryon production through the pro- 
cess e^e~ — >■ 7/^" — >■ Hqq/ (71) + Q + Q' versus the parameter 
Ptcut, where n stands for the intermediate diquark state. 



Ptcut 


GeV 


5 GeV 


10 GeV 


15 GeV 


20 GeV 




0.440 


0.420 


0.373 


0.310 


0.236 




0.214 


0.201 


0.174 


0.139 


0.101 


total [o-Ecc) 


0.654 


0.621 


0.547 


0.449 


0.337 


'^Hbc(PSi]3) 


0.607 


0.590 


0.546 


0.480 


0.397 




0.303 


0.295 


0.273 


0.240 


0.198 


'^Hf,J[iSo]a) 


0.441 


0.427 


0.390 


0.333 


0.265 




0.221 


0.214 


0.195 


0.167 


0.133 


total ((THi,^) 


1.57 


1.53 


1.40 


1.22 


0.99 




0.0486 


0.0468 


0.0423 


0.0355 


0.0273 


'^2,,{[lSol6) 


0.0233 


0.0223 


0.0196 


0.0158 


0.0116 


total (as, J 


0.0719 


0.0691 


0.0619 


0.0513 


0.0389 
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FIG. 3. The pt~ and y— distributions of the production 
channel e"*" + — !> 7 — >■ Stidn) + c + 5 at peak {y/S = 
mz), where n stands for the corresponding intermediate (6c)- 
diquark state. 



FIG. 4. The pt~ and j/— distributions of the production 
channel e+ + — s> — >■ H(,c(?^) + c + fe at peak 
= mz), where n stands for the corresponding interme- 
diate (6c)-diquark state. 



Total cross sections for the ^qq' baryon production at 
the Z'^ peak {^/S = mz) under several typical transverse 
momentum cuts, ptcut = 0, 5, 10, 15 and 20 GeV, are 
put in Table Hill With the increment of ptcut, the total 
cross section o'(Scc) decreases by 5%, 12%, 18% and 25% 
in a step-by-step way; the total cross section a (She) de- 
creases by 3%, 8%, 13% and 19% respectively; and the 
total cross section a{Ebb) decreases by 4%, 10%, 17% and 
24% respectively. It shows 5cc and Ebb have similar pt 
distributions, while Ebc is slightly different. 

Total cross sections for the Eqqi baryon production at 
the Z" peak (^5 ~ mz) under several typical rapidity 
cuts, Ucut = 2.0, 1.5, 1.0 and 0.5, are put in Table HVl 
With the decrement of Ucut, total cross section o'(Scc) 
decreases by 8%, 20% and 45% respectively; total cross 
section a{Ebc) decreases by 8%, 20% and 45% respec- 
tively; total cross section <7{Ebb) decreases by 5%, 18% 
and 43% respectively. It shows all the baryons have sim- 
ilar rapidity distributions. 

As a summary, 

• The SELEX collaboration [l^ has observed the 
S+ events, however, neither BABAR nor Belle 
collaboration has found the evidence for Sec in 
the related processes [H, H^. The CERN Large 
Hadronic Collider (LHC) or the future high lu- 
minosity e^e~ colliders (the super-Z factory, the 



TABLE IV. Total cross sections (in unit: pb) at the peak 
(\/S = mz) for the Hqq/ baryon production through the pro- 
cess e'^e~ — >■ 7/^" — >■ Sqq/ (n) + Q + Q' versus the parameter 
Vcut, where n stands for the intermediate diquark state. 



Vcut 


2.0 


1.5 


1.0 


0.5 


c^H^dSSila) 


0.424 


0.389 


0.310 


0.171 


'^E„(IiSole) 


0.206 


0.189 


0.151 


0.083 


total ((Th^^) 


0.630 


0.578 


0.461 


0.254 


'^Ht,([3Sils) 


0.593 


0.546 


0.434 


0.239 


'^Hbc(PSil6) 


0.297 


0.273 


0.217 


0.119 


'^3bc([lSol3) 


0.433 


0.400 


0.319 


0.175 


'^E,,([lSole) 


0.217 


0.200 


0.159 


0.088 


total ((TH(,^) 


1.54 


1.42 


1.13 


0.62 




0.0485 


0.0461 


0.0377 


0.0214 


0"S,,(IlSol6) 


0.0233 


0.0224 


0.0185 


0.0106 


total ((THoJ 


0.0718 


0.0685 


0.0562 


0.0320 



GigaZ project of ILC, and etc.) shall provide bet- 
ter platforms for measuring these baryons. Sev- 
eral studies on the heavy ba ryon pr oduction can 
be found in Refs. 0, i, 0, [H, iMll . Recently, it 
shows that even the triply heavy baryons could be 
measured at the LHC [33. 



5 



Fi/ra = 0.36 ±0.21. 
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Sizable events can be produced, which can provide 
a good platform for precise studies. More explic- 
itly, at the super- Z factory, we shall have 6.5 x lO'' 
5cc, 1.6 X 10^ She and 7.2 x 10'^ Sf,fc events per year. 

• To be a useful reference, several typical pt and y 
cuts have been adopted to show the baryon produc- 
tion properties in detail. These show that all the 
baryons have similar rapidity distributions, and the 
Sec and Sbh baryons have similar pt distributions, 
while the pt distribution of Sbc is slightly different. 
Even applying a large pt cut up to 20 GeV, wc can 
still have sizable events that are about half of the 
case without any pt cut. 

• As a final remark, there are two dominant decay 
channels S+ -> A+ K~ 7r+ and S+ p K~ 
for S+ . Setting Ti,2 to be the decay widths of these 
two channels, we have llSl: 



If taking ri/r2 = 0.36 and setting the relative pos- 
sibility for various light quarks as w : d : s ~ 1 : 1 : 



FIG. 5. The pt~ and j/— distributions of the production 
channel e+ + — >• 7/^" — >■ H6i,(n) + b + b at Z'^ peak 
(s/S = mz), where n stands for the corresponding interme- 
diate (fe6)-diquark state. 



• We have made a detailed study on the produc- 
tion properties of the doubly heavy baryon Sqq' 
through e'^e" annihilation. At the high luminosity 
supcr-Z factory and the GigaZ project at the ILC, 



, it IS tound that at the super- 
the luminosity C oc 10'^*cm~^s~^, we can obtain 7.4 
X 10'^ S+ events from the first channel and 2.1 x 
10^ S+ events from the second channel. 

Acknowledgements: This work was supported in 
part by the Fundamental Research Funds for the Cen- 
tral Universities under Grant No.CQDXWL-2012-Z002, 
by Natural Science Foundation of China under Grant 
No.11075225, No.ll275280 and No.ll205255, and by the 
Program for New Century Excellent Talents in University 
under Grant NO.NCET-10-0882. 



[1] G. T. Bodwin, E. Braaten, and G. P. Lepage, Phys. Rev. 

D51, 1125 (1995); 55, 5853 (1997). 
[2] J. Jiang, X. G. Wu, Q. L. Liao, X. C. Zheng, Z. Y. Fang, 

Phys. Rev. D86, 054021 (2012). 
[3] J. P. Ma and Z.X. Zhang, Sci. China, Pliys. Mech. Astron. 

53, 1947 (2010). 
[4] J. Erler, S. Heinemeyer, W. HoUik, G. Weiglein, P.M. 

Zerwas, Phys. Lett. B486, 125 (2000). 
[5] A. PetrelU, M. Cacciari, M. Greco, F. Mahoni, and M. 

L. Mangano, Nucl. Phys. B514, 245 (1998). 
[6] A. V. Berezhnoy, V. V. Kiselev, A. K. Likhoded, and A. 

I. Onishchenko, Phys. Rev. D57, 4385 (1998). 
[7] A. V. Berezhnoy, V. V. Kiselev, and A. K. Likhoded, 

Phys. At. Nucl. 59, 870 (1996). 
[8] D.A. Giinter and V.A. Saleev, Phys. Rev. D64, 

034006(2001). 

[9] M.J. Savage and M.B. Wise, Phys. Lett. B248, 177 
(1990). 

[10] N. BrambiUa, A. Vairo and T. Rosch, Pliys. Rev. D72, 
034021 (2005). 



[11] S. Fleming and T. Mehen, Phys. Rev. D73, 034502 
(2006). 

[12] T. Hahn and M. Perez- Victoria, Comput. Phys. Com- 

mun. 118, 153 (1999). 
[13] R. Kleiss and W. J. Stirhng, Comput. Phys. Commun. 

40, 359 (1986). 
[14] G. P. Lepage, J. Comput. Phys. 27, 192 (1978). 
[15] C. H. Chang, J. X. Wang, and X. G. Wu, Comput. Phys. 

Commun. 177, 467 (2007); Comput. Phys. Commun. 

181, 1144 (2010); X. Y. Wang, X. G. Wu, Comput. Phys. 

Commun. 184, 1070 (2013). 
[16] C. H. Chang, C. Driouich, P. Eerola, and X. G. Wu, 

Comput. Phys. Commun. 159, 192 (2004); C. H. Chang, 

J. X. Wang and X. G. Wu, Comput. Phys. Commun. 174, 

241 (2006); Comput. Phys. Commun. 175, 624 (2006); X. 

Y. Wang and X. G. Wu, Comput. Phys. Commun. 183, 

442 (2012). 

[17] J. Beringer et al., Particle Data Group, Phys. Rev. D86, 
010001 (2012). 

[18] M. Mattson et al, SELEX Collaboration, Phys. Rev. 



6 



Lett. 89, 112001 (2002); M. A. Moinester et al, SE- 
LEX Collaboration, Czech. J. Phys. 53, B201 (2003); A. 
Ocherashvili et al., SELEX Collaboration, Phys. Lett. 
B628, 18 (2005). 
[19] B. Aubert et al, BABAR CoUaboration, Phys. Rev. D74, 
011103 (2006). 

[20] R. Chistov et al., Belle Collaboration, Phys. Rev. Lett. 

97, 162001 (2006). 
[21] A. F. Falk, M. Luke, M.J. Savage, and M.B. Wise, Phys. 

Rev. D49, 555 (1994). 
[22] V. V. Kiselev, A. K. Likhoded, and M.V. Shevlyagin, 

Phys. Lett. B332, 411 (1994). 
[23] M. A. Doncheski, J. Steegborn, and M. L. Stong, Phys. 

Rev. D53, 1247 (1996). 
[24] S.P. Baranov, Phys. Rev. D54, 3228 (1996). 
[25] V. V. Kiselev and A. K. Likhoded, Phys. Usp. 45, 455 



(2002). 

[26] C. H. Chang, C. F. Qiao, J. X. Wang, and X. G. Wu, 

Phys. Rev. D73, 094022 (2006). 
[27] C. H. Chang, J. P. Ma, C. F. Qiao, and X. G. Wu, J. 

Phys. G34, 845 (2007). 
[28] J. W. Zhang, X. G. Wu, T. Zhong, Y. Yu, and Z.Y. Fang, 

Phys. Rev. D83, 034026 (2011). 
[29] J. P. Ma and Z. G. Si, Phys. Lett. B568, 135 (2003). 
[30] Z. J. Yang and T. Yao, Chin. Phys. Lett. 24, 3378 (2007). 
[31] S. Y. Li, Z. G. Si, and Z.J. Yang, Phys. Lett. B648, 284 

(2007). 

[32] Y.Q. Chen and S.Z. Wu, J. High Energy Phys. 08, 144 
(2011). 

[33] T. Sjostrand, S. Mrenna, and P. Skands, J. High Energy 
Phys. 05, 026 (2006). 



